Experimental observation of a complex periodic window 
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The existence of a special periodic window in the two-dimensional parameter space of an experi- 
mental Chua's circuit is reported. One of the main reasons that makes such a window special is that 
the observation of one implies that other similar periodic windows must exist for other parameter 
values. However, such a window has never been experimentally observed, since its size in parameter 
space decreases exponentially with the period of the periodic attractor. This property imposes clear 
limitations for its experimental detection. 
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The emergence of regular behavior is one of the most 
studied topics in nonlinear dynamical systems. It is 
known that by the changing of an accessible parameter, 
chaos pj and periodic 0] behaviors will be observed. 

The expectation of finding stable periodic behavior in- 
side chaotic regions in parameter space depends on the 
sizes and shapes of the parameter regions, regarded as 
periodic windows (PWs), for which stable periodic orbits 
(POs) are found. A PW is a region in parameter space 
that indicates parameter values for which one finds the 
lowest periodic attractor of period P, plus the period- 
doubling cascade with attractors of period P 2n , with 
n EN (|. 

For systems whose chaotic attractors have only one 
positive Lyapunov exponent as the Chua's circuit, con- 
sidered in this experiment, a special type of PW, regarded 
as complex periodic windows (CPWs), is everywhere ob- 
served in parameter space. The appearance of one such 
window implies in the appearance of an infinite num- 
ber of self-similar others that appear side by side aligned 
along a direction. In addition, CPWs have an extended 
characteristic in the parameter space. They visit large 
portions of the parameter space; i.e. one can still stay in 
the same periodic windows even if especially large varia- 
tions in two control parameters are made. Due to these 
two characteristics an arbitrary change in only one acces- 
sible parameter can replace chaos by periodic behavior, 
or vice versa. So a better understanding of a CPW is 
relevant to applications that relay either on a robust pe- 
riodic oscillation, as mechanical machines, or on a robust 
chaotic system, as chaos-based communication 

These CPWs, regarded as shrimps were extensively 
studied in maps jo, 0] and in periodically forced maps 
@,@| • However, only recently were these windows numer- 
ically observed in systems of ordinary differential equa- 
tions [l(J EH ■ The reason is that the parameter interval 
length, AV, of a CPW scales exponentially with — P, 
where P is the period P of the lowest-period periodic at- 
tractor of the CPW @| . Since CPWs have usually higher 
P, they are too tiny to be observed, even though these 
tiny windows are extended in parameter space. 



This exponential scaling clearly imposes limitations on 
the experimental detection of such a periodic window, 
and arguably due to that, they have never been exper- 
imentally reported. However, for the Chua's circuit, it 
was numerically shown in Ref. [ll[ that such CPWs pos- 
sessing a low value for the lowest-period periodic attrac- 
tor (P=4) exist. This work is dedicated to experimen- 
tally report, for the first time, such a CPW. 

To certify that we observed a CPW, we show that 
there exists curves in parameter space where the POs 
are super-stable, and that these curves cross transver- 
sally at least twice, a necessary condition that defines a 
CPW. These parameter curves are detected by the indi- 
rect method of noting the parameter values at which the 
symbolic sequences, encoding the type of POs existing 
within the CPW, change. 

The well known Chua's circuit is shown in Fig. QJA). 
The control parameters are Ri = Rio — ARi and i?2 = 
R 2 q — AR 2 , where Rio and R20 have fixed values, ARi 
and Ai?2 are varied by precision potentiometers, with 
steps of 50mfl and 200mf2, in the ranges [0, 17]f2 and [1, 
5.5]fi, respectively. We obtained time series by recording 
the Vc\(t) voltage with a 12-bit analog-to-digital con- 
verter (ADC) at the rate of 400 ksamples/s. All the 
attractors were reconstructed by the Takens method [12] 
with time-delay t = 45.0 [is, which corresponds to 18 
data points. Then, the reconstructed attractors are made 
discrete by measuring Vci(t + r) when the reconstructed 
trajectory reaches the section Vci(i)=~2.25V in a clock- 
wise orientation. The value of Vcxit+r) when the recon- 
structed trajectory realizes its nth crossing in this section 
is denoted by V^ 1 . 

In Fig. QjB), we show the parameter space of this 
circuit. There, solid black circles represent parameter 
values for which one obtains the lowest period PO. Along 
the left border between chaos and the PW [parameters 
indicated by letters "a" within the boxes of FigQjB)] in 
these two PWs, chaos is replaced by a stable (period-3 or 
period-4) attractor by a tangent bifurcation by increasing 
Ai?2- In the other borders, [parameters indicated by 
letters "b", "c" and "d"], the lowest-period PO inside 
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Figure 1: (A) Scheme of Chua's circuit. Their component val- 
ues are: R w » lAkQ, R20 « 37ft, Ci « 4.7nF, C 2 » 56nF, 
L « 9.2mH. (B) Parameter space for the experimental Chua's 
circuit showing the period-3 and period-4 windows. Filled 
black circles represent parameters for which the lowest pe- 
riod POs are observed, filled light gray circles represent the 
higher period PO that appear by periodic-doubling bifurca- 
tions, and filled dark gray circles represent parameters for 
the closest chaotic attractor to the PWs. The straight lines 
indicate parameter values for the sets of time series X and Y. 



the PWs bifurcate and chaos (outside the PW) is reached 
after a period-doubling cascade by modifying AR2 . 

To illustrate our analysis techniques, we first use a 
symbolic representation to characterize the lower period 
POs that appear for the parameters nearby the borders 
between the PWs and the chaotic regions, indicated by 
the letters a, b, c and d, in the boxes I, II, and III, in 
FigQjB). We use data sets collected varying Ai? 2 along 
the lines X and Y for ARx= 3.0ft and AR 1 = 12.5n, 
respectively, as shown in Fig.QjB). 

The symbolic characterization of these POs is done by 
encoding them by the approach in Ref. [13], using the 
properties of the nearby chaotic attractors. The return 
maps of the reconstructed chaotic attractors for param- 
eters in the borders a, b, c, and d in box III, are shown 
in Figs. [2]JA-D). These maps as well as the other chaotic 
attractors at the borders in both period-3 and period-4 
windows display return maps typical of either unimodal 
(one maximum) or bi-modal (one maximum and one min- 
imum) maps, and they can be partitioned by the critical 



points. The partitions are in the maximal and minimal 
points, assigned by Vx and Vi. So a trajectory point in 
the interval Vex < Vx is encoded by '0', a trajectory point 
in the interval V\ < Vex < V2 is encoded by 'V, and a 
point in the interval Vex > Vi is encoded by '2'. A stable 
period- P orbit can be encoded by comparing its mapping 
with the mapping of the nearby chaotic attractor, and de- 
pending on the position of the POs points with respect 
to the partition points, a PO can be encoded by a se- 
quence &\a% ■ ■ ■ sp, where is a symbol of the alphabet 
Si = {0, 1, 2}. For chaotic attractors close to the borders 
with the period-3 window, in box I, the chaotic returning 
maps are uni-modal, with only one critical point V±. In 
the side a of the window, in box I, we obtain the symbolic 
sequence 101 and in the right side b, the sequence 100. 
All the POs in the left side of this window are encoded by 
101 and the ones on the right side by 100. The period-4 
POs, in the period-4 window, close to the borders a, b, 
c and d, in box III, [whose return maps can be seen in 
Figs. [U[A-D), respectively] are encoded by the sequences 
1001, 1000, 2000 and 2000, respectively. 

In fact, as one varies a control parameter, the symbolic 
sequence of a stable PO changes if some periodic point 
crosses a critical point of the return map hj|. This mech- 
anism is responsible for the changes in the symbolic se- 
quences of the stable POs in the period-3 window. There, 
the symbolic sequence 101 changes to 100 when the PO 
crosses the critical point V% . 

We name £ the return map of the closest chaotic attrac- 
tor to the period-P PO, and O a stable PO with points 
Vq , ■ ■ ■ , Vq . Assuming that the return map £ can be 
used as an approximation to calculate the first derivative 
of the orbit points of a PO inside a PW, then the orbit 
O is stable if 

A < 1 (1) 

with A = I YiiLi d0i~\- If a contains a critical point, 

a point on the extremum of the map, A = 0, and we 
say such an orbit is superstable. For parameters e-close 
to a parameter for which a super-stable PO exists, Eq. 
(U) is satisfied, which means that it exists a PW in the 
neighborhood of parameter lines for which V c = V\ . 

A similar mechanism governs the changes in the sym- 
bolic sequences of the stable POs inside the period-4 re- 
gion. The difference now is that we have two critical 
points, Vx and V2, which makes Eq. ([1]) satisfied in pa- 
rameter curves for which either V c = Vx (which defines 
the critical curve Svx) or Vq = V% (which defines the 
critical curve SV2), or V c = Vx and V c = V%. It is 
typical for this type of CPW that the PW appears not 
only for the parameter point for which V c = Vx and 
Vq — V2, a zero measure point in parameter space, but 
also along the curves Svx or Sv2- These two curves form 
the spines introduced in Refs. 0, Q. 

Three important characteristics grant to this window 
the status of being a CPW: (i) if there is one CPW, then 
a countable infinite number of others must exist, with 
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Figure 2: [Color online] Return maps (black points) of the 
Poincare section of the chaotic attractors obtained using the 
parameters indicated by the borders a, b, c and d in box 
III of Fig. [TJ respectively Figs. (A)-(D). We also show the 
return maps (blue circles) of the periodic attractors obtained 
for the closest parameters to these borders. The vertical lines, 
passing through the maximum and the minimum define the 
partition points. In (D') is shown a zoomed view of the min- 
imum of the return map in (D). 
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Figure 3: [Color online] (A) The encoding of all the period-4 
POs found in the CPW. (B) Sketch of the critical lines (Svi 
and SV2) structure of the CPW, disregarding the existence of 
characteristic (iii) that causes the appearance of structures as 
illustrated in Fig. [4] 



sizes that decreases exponentially [Eq. J5])] as the period 
of the POs increase, (ii) The two critical curves Svi and 
Sv2 cross transversally at least twice. For the parameters 
where the crossings happen, the PO has an orbit point 
Vc x = V\ and another Vc\ = V%. (iii) POs with the same 
period coexist. 

Concerning characteristic (i), for quadratic maps one 
should expect that 

AV{P) oc e~ pp (2) 

as shown in Ref . @ , with V being the parameter interval 
length of a CPW, and P the period of the lowest-period 
periodic attractor. Also, from Q, we have that f3 = 2Ht, 
where Ht is the topological entropy or Lyapunov expo- 
nent of the bordering chaotic region 0, But, in fact, 
for flows such as the Chua's circuit containing Shilnikov's 
homoclinic orbits [l5| , a two-parameter analysis per- 
formed in the neighborhood of this orbit shows that it 
exists a countable (infinity) number of CPWs that ap- 
pear side-by-side in parameter space following the same 
exponential scaling law that describe the appearance of 




Figure 4: Illustration of the structures that might appear in 
a CPW due to characteristic (iii). 



the homoclinic orbits. This exponential scaling law is of 
the form of Eq. (J2j) , and as shown in Ref. [13], /? = tt^, 
with p and u> representing the real and imaginary parts of 
the eigenvalues of the focus point associated with the ho- 
moclinic orbit responsible for the generation of the many 
CPWs. 

We estimate that for this experimental circuit f3 ~ 2, 
in Eq. ^ , for a parameter region in the vicinity of the 
observed period-4 CPW. That means that in order to 
observe a higher period CPW, with period Ph = 4 2n , 
with n G A/", associated with the observed period-4 CPW, 
we should have a potentiometer with a resolution (step 
size) of 8Apexp(~ 2 ( p,l ~ 4 )\ 8 being roughly an average 
width of the period-4 CPW observed. So, in order to 
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observe a period-8 CPW, we would need a potentiometer 
with a resolution of about 0.14mf2, which is much smaller 
than our experimental resolution. Numerical simulations 
realized in a similar Chua's circuit, reported in Ref. flll | . 
show that CPWs with attractors of period lower than 4 
exist. However, their sizes are of the order of 20 times 
smaller than a period-4 large CPW, similar to the one 
observed experimentally. Therefore, for the resolution of 
our experiment, we do not expect to find the many others 
numerically found CPWs, but only this "giant" one. 

To detect the existence of the critical curves, we search 
for transitions in the symbolic sequence of the POs closer 
to the borders between the PW and chaos. In box II, the 
PO encoded by 1001 at the border a changes its encoding 
to 2001 at the border b. So, between these two borders, 
there is a parameter Ai?2 for which at least one point 
of the period-4 orbit is — Vi. Thus, within these 
borders there must exist a curve Sv2- In box III, the PO 
encoded by 1001 (border a) changes its encoding to 1000 
(border b), indicating that within these borders there 
is a PO that visits the critical point V\. Thus, within 
these borders there must exist a curve Svi- In box III, 
the POs in both borders c and d are encoded by the 
symbolic sequence 2000, which suggests that within these 
two borders there must exist either (or both) curves Svi 
or Sv2- 

As we go from one side of the CPW to the other side by 
changing Ai?2 , for a fixed ARi , the points of the return 
map of the POs wander along an imaginary smooth curve 
This imaginary curve changes its form smoothly, as 
we vary Ai?2- For a AR2 close to a parameter where 
chaos is found (close to the borders a, b, c or d), £' re- 
sembles the return maps £ of the chaotic attractors. The 
curve £' can be constructed using all the POs observed 
in this CPW, for a constant Ai?i. Then, we estimate the 
location of the critical points of which provide us the 
encoding for the period-4 POs within the CPW, in Fig. 
[3](A). The curves Svi and Sv2 are located where two 
different colors (which describe the different encodings) 
meet. A curve Svi is the border line between two re- 
gions representing different encodings. Either '1001' and 



'1000', or '2001' and '2000'. A curve S V 2 is the border 
line between the regions that encode either '1001' and 
'2001', or '1000' and '2000'. Note that these curves cross 
transversally at least twice inside the windows, at the 
points where the regions that encode the four different 
types of POs meet. This is characteristic (ii) of a CPW 
Q. It can be understood by the way CPWs appear in the 
parameter space. The process can be described as hav- 
ing a normal PW which contains two curves Svi and Sy% 
that do not cross. One can imagine that both curves have 
a parabolic shape appearing side by side. As one changes 
a parameter of the circuit, the curve SV2 approaches Svi 
crossing it in at least two points forming a structure simi- 
lar to the one shown in Fig. [3KB), a sketch of a simplified 
version of what it could be really happening inside the 
CPW. There, one sees that some regions in the param- 
eter space that represent POs with some encoding (e.g. 
'1001') do not border a region with some other encoding 
('2000'), except for the point where the curves Syi and 
Sv2 cross. And when that happens (excluding the atyp- 
ical case when the curves are tangent), there has to be 
at least one more crossing inside the CPW, so that the 
POs appear side by side other allowed POs. The rule is 
'1001' appears aside '1000', which appears aside '2000', 
which appears aside '2001', which appears aside '1001'. 

Such a rule can be apparently violated due to char- 
acteristic (iii) which leads to points where two or three 
different regions meet, as represented in Fig. [H But note 
that, in fact, the line Svi does not cross the line Sv2, and 
thus, the rule that describes the crossing between these 
lines is not violated. Internal noise and parameter fluctu- 
ations of the circuit partially destroy the CPW. Adding 
the fact that we have limitations in our parameter resolu- 
tion, we do not expect to identify all these fine details of 
the CPW, but rather a lower resolution picture, in which 
this rule might be apparently violated. 
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